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Abstract 



Reviewing the semiclassical theory for the parametric level density fluctu- 
ations, we show that for large parametric changes the density correlation 
function, after rescaling, becomes universal and coincides with the leading 
asymptotic term obtained from Random Matrix Theory. The advantage of 
the semiclassical approach is to provide a simple recipe for the calculation of 
the non-universal scaling parameter from elements of the underlying classi- 
cal dynamics specific to the system. We discuss recent improvements of the 
theory, which introduce a requantization of the smoothed level density. 

PACS numbers: 05.45. +b 
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I. MOTIVATION 



One of the interesting new developments in the study of quantum manifestations of 
classical chaos concerns parametric correlations. For over a decade, most of the works on 
quantum chaos concentrated in spectral fluctuations of a fixed Hamiltonian. This large effort 
led to a solid body of evidence that if a system is chaotic in the classical limit, its quantum 
spectrum exhibits universal fluctuations which coincide with those of a suitable ensemble 
of random matrices M. There was also some parallel activity in trying to understand the 
statistical behavior of quantities related to the response of a chaotic Hamiltonian to external 
perturbations, such as Landau-Zener transitions [TJJ and energy level curvature ||. But it 
was only a few years ago that Szafer, Simons, and Altshuler ||,|5] showed that a system 
whose spectrum follows closely the universal fluctuations predicted by random matrix theory 
(RMT) will also present universal parametric behavior. More precisely, these authors have 
concluded that when a chaotic Hamiltonian depends on some external parameter X, any 
correlation function of spectral fluctuations taken at different values of X becomes system- 
independent after a proper rescaling. 

In this paper we study some universal spectral features of a quantum mechanical prob- 
lem whose Hamiltonian H generates chaotic motion in its classical limit. Considering a 
parametric dependence of H on X, the Schrodinger equation reads 



In what follows, the parameter X will be associated with any tunable physical quantity, 
provided that its variation neither alters significantly the classical dynamics of the system, 
nor breaks or restores any symmetry. The first condition is easily realized in the semiclas- 
sical limit, since classically small parameter variations are responsible for significant level 
fluctuations. The second condition is not difficult to be satisfied in most cases of practical 
relevance. 

Previous works have concentrated on the level velocity correlation function 



where A denotes the mean level spacing around the state v. This function, although straight- 
forward to evaluate numerically, is not suitable for any exact or systematic analytical evalu- 
ation. We would like to stress that the difficulty is intrinsic to the function, since it requires 
the precise knowledge of every level as a function of X. Therefore, instead of the usual 
2-point Green's function, one would need a combination of arbitrarily high iV-point Green's 
functions, which is very impractical. A few years ago, Berry and Keating J7| proposed a 
semiclassical approximation based on the density correlation function for C(X). Their result 
agrees very nicely with the numerical results [§] for the large X limit, indicating that the 
2-point Green's function gives a good asymptotic expansion. For the small X region, their 
analytical expression gives, as one would expect, very poor results. 

For these reasons, in order to study parametric statistics semiclassically, we will rather 
focus on a formally more amenable quantity, namely, the density correlation function 




(1) 





2 



E,X 

where the fluctuating density p? 1 = p — (p) is defined in terms of the level density 

and the average density (p), or inverse mean level spacing 1/A. Sometimes, in the semi- 
classical approach, it is convenient to use the cumulative level density, N(E) defined as 
N(E) = Jf oc dE'p(E'), instead of the density itself as we shall see in Section IV. 

All statistical measures of spectral fluctuations presuppose a constant average level 
spacing, which is not guaranteed for general families of Hamiltonians. Indeed, the mean 
level spacing for smooth Hamiltonians is accurately approximated by the Weyl rule: 
N Weyl (E, X) m V(E,X)/(2nh) d , where d is the number of degrees of freedom and 

V(E,X) = Jd d pd d qe(E-H(p,q,X)) , (5) 

G being the unit step function. It is therefore desirable to work with the "volume spectrum" 
rather than directly with the energy spectrum, as proposed by Goldberg and collaborators 
0j. So as to avoid some technicalities in the classical theory, we will adopt the alternative 
procedure of defining an "unfolded classical Hamiltonian" 

H'(p,q,X) = v(H(p,q,X),X) . (6) 

The classical motion is then identical to that of the original Hamiltonian, except for a 
rescaling in time. Quantizing (||) we obtain levels E' v that fluctuate with X about stable 
mean values u(27ih) d , i.e., for the unfolded system, A = (2nTi) d . (The primes will be omitted 
in what follows.) 

The paper is organized as follows: In Section II, we present a reminder of the main 
parametric random matrix results, extracting the relevant asymptotic correlation functions. 
In Section III the standard semiclassic theory is discussed, putting emphasis on the diagonal 
approximation and its sequels. In Section IV we discuss the implications for K (fl, X) of 
the recent proposal of Bogomolny and Keating [|Kj for requantization of the smoothed level 



density. We summarize the results and conclude with a general discussion about the range 
of validity of the semiclassical approach in Section V. 



II. REMINDER OF THE RANDOM MATRIX RESULT 

The random matrix calculations leading to an exact expression for K (Q, X) have become 
relatively standard. Several references can serve as introduction to the subject (see, for 
instance, [pTT|— [13|1 ) . The detailed derivation of K(fl,X) was first given in Ref. |J in the 



context of disordered metallic systems. The technique employed for such continuous systems 
can be easily adapted to discrete Hamiltonians. It consists in expressing the one-point 
Green's function Gji(E) = (l\(E — H+iO + )~ 1 \j), where / and j label arbitrary basis states, in 
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terms of an integral over N commuting and N anticommuting variables, where N represents 
the number of basis states, and then directly performing the average over a certain Gaussian 
ensemble of Hamiltonians, {H} NxN . The quartic term thus generated can be decoupled 
through a Hubbard-Stratonovich transformation at the expense of introducing a graded 
matrix, also known as a supermatrix. The size and symmetry of this supermatrix depends 
on the symmetry class of the Hamiltonian and the order of the correlation function in 
question. For instance, for the Gaussian unitary ensemble (GUE) the calculation of the 
density correlation function requires an 8 x 8 supermatrix transforming according to the 
group U(l,l|2). After integrating out the commuting and ant i- commuting variables, one then 
has to carry out an integration of an effective action that involves the supermatrix elements. 
This is usually done in the saddle-point approximation, exploring the condition N ^> 1. 
Notice that this technique greatly simplifies the problem: The 2N variables are reduced 
to a small, workable set. However, we should remark that the apparently straightforward 
steps just described involve rather difficult mathematical considerations. For the important 
details the reader should consult the references cited above. 

One of the most important conclusions drawn from Refs. []5],|6| is the existence of a single 
quantity controlling the scale of parametric fluctuations of all spectral functions of chaotic 
Hamiltonians. This result is obtained explicitly in the supersymmetric formulation of the 
disordered metallic system in the diffusive regime (the so-called zero mode), as well as (but 
less surprisingly) in RMT when N — > oo. More precisely, if X is the external parameter, 
the rescaling x = X/X c , where 



X r 



1 dE v (X) 
A dX 



n2v -i/2 



(7) 



eliminates all system dependence in the parametric correlation functions, up to a global 
prefactor. We will later argue that a similar conclusion can be obtained from a purely 
semiclassical approach. 

We now present the results of Ref. || in the rescaled form, assuming also that all energies 
are expressed in terms of A, introducing uj = Q/A. In the absence of time-reversal symmetry, 
the final result for the density correlation function is 



I rl POO 

- Re / d\ / d\\ exp 
2 J-i Ji 



k GUE (uj,x) 

where k = Kj (p) 2 and the free energy is given by 



—F 



GUE 



(A, 



^GUE 



7l 2 X 2 



(A?-A 2 )+mcu+(A-A 1 ) 



(8) 



(9) 



with uo + = uj + i0 + . The integration in (§) can be carried out explicitly, and gives 



k 



— - Im I erfc f - 
Ax 2 \ \ 



V2 



erfc 



uj /x + mx \ / — uj jx — inx 



V2 



s/2 



(10) 



where erfc is the complementary error function [ I4"| . 
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For systems where time-reversal symmetry is preserved, modelled by the Gaussian Or- 
thogonal Ensemble (GOE), the density correlation function is 



i roc roc (1 - A 2 ) (A - Ai A 2 ) 2 exp —F GUE (X, Ai, A 



k a0E {u,x) — Re Jd\ ^ ,/V — _^ _ _ _ _ - , 

with the free energy 

f goe = ^(2A 2 A 2 - X\ - A 2 - X 2 + 1) + iiru+{X - X 1 X 2 ) . (12) 

In distinction to the previous case, here the triple integral in ([11]) cannot be expressed in a 
closed form. 

We will now focus on a regime when one either has x ^> 1 or u ^> 1. Considering initially 
large variations of the external parameter, one finds that the leading contributions to the 
multiple integrals come from the regions around two points: A = Ai = A2 = 1 (A) and 
A = — Ai = — A2 = —1 (B). In passing, we remark that there exists a mapping of RMT into 
the problem of many particles interacting by a 1/r 2 potential in one dimension |15| . In that 
context, the contribution (A) gives the "hydrodynamical" , q = limit of the particle density 
correlation function, whereas (B) is related to the "Friedel", q = 2kp oscillations. Here, (A) 
is related to the monotonic, non-oscillating decay of the correlation function; (B), on the 
other hand, yields a series of oscillating terms in to. To obtain an asymptotic expression 
for k(u,x), we expand the free energies around these points and retain the (linear in A's) 
lowest-order terms. Performing the exponential integrals (see the Appendix [A] for details), 
we find that 

k GUE (u x) ~ -- ^ - — + - C ° s(27ra;) (13) 

^' X) 2 (TT 2 ^ 2 + 7T% 4 ) 2 + 2 7T 2 ^ 2 + 7T% 4 ^ 

and 

k GOE (u x) - ^ ~ + 1 cos(27ro;) 

K , ) (vr 2 cu 2 + vr 4 x 4 /4) 2 2 (vr 2 ^ 2 + tt 4 x 4 /4) 2 ' 1 } 

The right-hand sides of both equations represent the leading terms of an asymptotic expan- 
sion on both x and uo whenever \tt 2 x 2 j3 / 2 — muo\ ^> 1, with (3 = 1 for the GOE and (3 = 2 for 
the GUE. At x = one has the usual 1/uj 2 decay of the density fluctuations, plus the first 
oscillating terms. Notice also that the oscillating and non-oscillating terms are held only 
to their lowest order in x~ 2 and u~ 2 . For the GOE, this means that we are not showing 
a non-oscillating term of order u~ A . For the GUE, higher orders terms can be obtained 
without much effort by recalling (|T0|) . 

The energy oscillations present in Eqs. (|T3| ) and ( |T^D are related to the quantum (discrete) 
nature of the energy spectrum. For the GUE case they appear even at the lowest order in 
u>~ 2 ; in fact, for x — 0, the asymptotic form (|I3J) is actually exact and applies to all to scales. 
For the GOE case, since level repulsion is less pronounced, the oscillations only show up at 
the next-lowest order, namely, they are weaker corrections to the leading uo~ 2 decay. 

One is tempted to imagine that the uo oscillations should not be present in a formulation 
where there is no intrinsic energy level quantization. Indeed, in the following section we 
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will verify that the semiclassical trace formula (applied in the usual way and based on the 
assumption that the underlying classical dynamics is chaotic) only yields the monotonic 
decay of k(u,x). Still, a problem emerges as to whether one can systematically obtain 
oscillating terms, or any parametric correlation function beyond the asymptotic limit, in a 
semiclassical framework. This question was partially answered by Bogomolny and Keating 
flCfl when they calculated k(u, x — 0) using an expansion in periodic orbits and forced energy 
level quantization as an additional constraint. We will get back to this point later. 

Another important aspect, often neglected, but particularly relevant to mesoscopic phe- 
nomena, is the applicability range of the semiclassical approach used in the study of para- 
metric variations. For this purpose, the calculations below will be carried out with special 
care. 



III. SEMICLASSICAL APPROACH 

The fluctuating part of the level density p^ l {E, X) can be smoothed over an energy range 



rj/Ti and expressed, in the semiclassic limit, by the Gutzwiller trace formula |16 



P fl (E,X) = JL E , A ^ ( J: nil/2 exp (Us, - ^ - V -\r\T 7 ) , (15) 



27th %\det(M'-I)\V* *\h 7 2 7 ft 1 

where T 7 , M 7 , S" 7 , and v 1 stand, respectively, for the period, monodromy matrix, action, 
and Maslov index of a periodic orbit labelled by 7. The sum is performed over all primitive 
orbits 7 and their repetitions r (positive and negative), irrespective of multiplicity due to 
time-reversal or other symmetries. The smoothing of p* 1 can be done in other ways than the 
exponential cutoff introduced in ([15]). Later on, we shall also cutoff the sum over periodic 
orbits sharply at a certain period T*. Both h/T* and r]/2n are always taken larger or equal 
to A. 

Let us write the correlation function K (fi, X) in terms of the semiclassical level density 
as given by Eq. fll5D , 

K *cr n x) = —L-( y A r (E + -,X + -)jL, r ,(E--,X-- 



x exp <j — 



7,r,7',7" 

i 



2' 2 J 7 V 2' 2 J 



+ i^(ru 7 -r'v Y )- 1 l(\r\T 7 +\r / \T 7 ,)}j _, (16) 

where A ir = T 7 | det(M^ — J)| -1 / 2 . In order to further proceed analytically, we restrict the 
calculations to a regime where classical perturbation theory is applicable. This means that 
we shall consider variations of E and X which classically imply a small change in the actions 
S T Such a change can nonetheless be very large in the scale of h, corresponding to large 
quantum effects. Under this condition, it is a good approximation to write 

7 V 2 ' 2 J 7V ; dE 2 dX 2 

= S,(E,X)±T,(E,X)^±Q, / (E,X)^ + -- - , (17) 
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which defines Q 1 as the parametric velocity of the orbit action. The general expression 
for the parametric variation of the action is presented in Appendix B. In the absence of 
any dependence on Planck's constant, one can neglect the energy (or parameter) correc- 
tions to as compared to the strong energy dependence of the exponential term. In this 
approximation, consider the evaluation of 

Iw(T) = L^iE^-S^X)]^ (lg) 



For orbits with period shorter than the Heisenberg time = h/A, one does not in general 
expect to find pairs of orbits (7, 7') with actions differing by less than h, unless they are 
symmetry related orbits. As a consequence, upon energy averaging, J 77 /(T < tu) ~ ^7^7,7') 
where g 7 is the multiplicity of symmetry-related orbits. (We take here the simplest case 
where any symmetry holds for all considered parameter values.) This is the essence of the 
diagonal approximation. On the other hand, due to the exponential proliferation of orbits, 
for sufficiently long times, i.e., for a small smoothing parameter r], one enters a regime 
where there is an abundance of pairs of periodic orbits (7, 7') satisfying 5 7 — Sy < %. 
This situation is expected to occur for time intervals longer than tjj (or, correspondingly, 
for E < A). Under these circumstances, the diagonal approximation is no longer valid. 
Actually, for such long times even the validity of the trace formula is problematic, so that 



the unsmoothed energy spectrum should be resummed |Tj]. At this point the necessity of 
smoothing becomes evident: it keeps the approximations under control. 

Restricting ourselves to the range where the diagonal approximation is valid, we have 

Kb c (V,X) = — (j2g,\A ir (E,X)\ 2 



7,r 



X 



exp{~ \rT 1 (E,X)n + rQ^(E,X)x] -2~|r|T 7 }\ _. (19) 
I n 1 J n J / e,x 



To evaluate the average in ([19]), we first notice that, in a fully chaotic regime, the proportion 
of high period orbits that are repetitions of lower periods is exponentially small, so we neglect 
all \r I ^ 1. In the next step we recast fllTf ) as an integral by defining the smooth interpolating 
functions A 2 (E,X,t) through the relation 

f dtA 2 (E, X, t) £ 5{t - T 7 ) = £ A 2 (E, X) . (20) 

77 

Further simplification can be achieved using the uniformity principle over periodic orbits 
(also known as the Hannay-Ozorio de Almeida sum rule) 



A 2 (E,X,t)Y,Ht-T,) > \t\. (21) 

It is important to notice that this result independs on the "center-of-mass" variables E and 
X. We can then evaluate 

K£(Q, X) = ^pRe j™dt \t\ e*n/»-^l/» (e^W*) . (22) 
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The lower limit r is the period of the shortest periodic orbit. Of course it is simplistic to 
extrapolate the uniformity principle this far, but it is easy to improve the theory by including 
a few short orbits individually and hence increasing r. 

The average in Eq. fl22|) is evaluated over periodic orbits with fixed period t. In view of 
the fact that the rest of the integrand oscillates rapidly in time, it is important that the finite 
range in E and X allows many orbits with a given large period. There being two parameter 
families of periodic orbits, one should take Q^(E(X), X, t) and subsequently average only 
over X, or vice-versa. For the energy correlation with no variation of parameters, there 
will be isolated periodic orbits of given period in the energy range where the averaging is 
carried out. Thus, for fixed X, the average in Eq. (^) runs over a large number of periodic 
orbits with parametric velocity Q 7 (t), which may be considered as samples of the probability 
distribution P t {Q). The latter is assumed to be Gaussian, 



p t(Q) = i exp 

'27cQ 2 (t) 



I g 2 
2Q 2 (t) 



(23) 



because our theory involves arbitrary parametric variations H(X). Even if there exist cor- 
relations among the actions of the periodic orbits 7(t), as postulated by |19 , these will in 
general be broken by the parametric velocities 



r 1 -! OH / 
Ql = Jo ^5x(pW'^)' X 



(24) 



evaluated along each orbit. The Gaussian width is a function of X and E (the center of the 
energy range), being defined as 



, 1 



Q 2 (t) 



iv(t)^ g ' w ' 



(25) 



7(t) 



where N(t) is the number of periodic orbits with period t. Note that we take Q = 0, since the 
shape of the shell for the unfolded system is affected by X, but not its volume, so dH/dX 
averaged over the shell is zero. Although the Gaussian assumption is not quite at a par 
with the uniformity principle, the arguments presented above make it plausible. Nothing 
is yet known of the relation between higher moments Q n and Q 2 , so as to obtain a more 
rigorous result. Numerical simulations for the particular case where the external parameter 
is a magnetic flux line || provide further support to (]23|) . 
The average in Eq. (|22| ) can be evaluated using (p3|): 



iQX/K\ 



exp 



X 5 



2h- 



Q 2 (E,X,t) 



(26) 



As we commented previously, the energy average was already utilized in obtaining the pe- 
riodic orbits with fixed period. Strictly, we should now average (|^) over X, but both Q 2 
and the Gaussian (|23| ) are smooth functions of A, so we may just evaluate the latter at the 
center of the averaging range. As for the time dependence of the mean square parametric 
velocity, we refer to previous treatments P,PU[, leading to the result 
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roc I r)TJ f)J-f \ 

W=a\t\, with a(E,X) = 2j dt ( — (p(t),q(t)) — (p(0),q(0))j . (27) 

^ \ / e.s. 

Here, the average over the periodic orbits has been substituted by an average over the 
entire energy shell, in accordance with the uniformity principle. It is important to note that 
the decay of the classical correlation function in fl2"7| ) need only be integrable, since many 
available chaotic systems do not exhibit full exponential decay of the correlations. Strictly, 
the linear dependence of Q 2 on t should only hold for times longer than that of the decay 



of the correlation function, but numerical investigations @,^1|] have found this feature to 
be quite robust. (Notice that there is a factor 2 missing in the definition of a in Ref. |§. 
Modifying a accordingly the semiclassical prediction for X c changes by a factor y/2. As a 
result, the agreement between the numerics and the semiclassical theory becomes excelent.) 

Substituting fl2~?|) and fl2"6| ) into Eq. ( p2[ ) and evaluating the integral in the limit r — > 0, 
we finally arrive at the main result of this section 

[9? + (aX 2 /2h + 2r]) 2 

where j3 = 2/g. The nongeneric long-wave oscillations in the correlation, characteristic of 
each individual system, appear as we choose r equal or greater to the period of the shortest 
periodic orbit. The above expression corresponds to the leading non-oscillatory term of the 
asymptotic expansion of the correlation function exactly calculated within RMT [Eqs. ( |T3D 



and and hence presumed to hold for generic chaotic systems. 

On the other hand, this formalism gives a recipe for calculating a, which is a system- 
specific quantity and fixes the value of the parametric rescaling factor X c . From this point- 
of-view, the semiclassical method and the RMT are complementary for a full understanding 
of K(Q,X). An important immediate application concerns possible transitions from GOE 
families to GUE families of Hamiltonians. If this occurs for a small variation of a second 
parameter as in |2U||, the classical parameter a should remain unchanged, whereas X c must 



grow because of increased level repulsion. By comparing the invariant form of Eq. ( |28|) with 
(|l3|) and fll4f) , we can predict that X c will increase by a factor \f2. 



IV. REQUANTIZATION AND ITS IMPLICATIONS FOR THE DENSITY 

CORRELATION FUNCTION 



In a recent paper, Bogomolny and Keating |L0| proposed a new semiclassical procedure 
to obtain the density correlation function K(Q, 0) and its parametric extension, K(Q,X). 
Their result for K (Q, 0) coincides with the exact GUE expression for systems with broken 
time-reversal invariance and gives the leading GOE oscillatory and non-oscillatory terms 
in powers of u~ 2 for time-reversal symmetric systems. We will here rederive their results 
within the wider context of parametric variations. 

The starting point of the derivation in Ref. Ill is standard and considers the sum over 



periodic orbits in the Gutzwiller trace formula (^) up to some value of T*. As a consequence, 
the information about energy scales smaller than h/T* is washed out, and strictly speaking, 
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the connection to RMT results cannot be made in such an energy range. However, one 
can still introduce information related to the discrete nature of the spectrum through the 



following "requantization" condition p2 



NT.(KVrj)=v + ±, (29) 

by which, for a given value of T*, one has a well-defined procedure to obtain a set of 
eigenvalues {E U (T*)}. It has to be stressed that such a procedure does not imply any 
control on accuracy. This is equivalent to say that, by increasing the value of T*, it is not 
guaranteed that {E U (T*)} will converge to the exact spectrum. Actually, this quantization 
rule is not applicable for T* > tu since for such periods the cumulative semiclassical level 
density can become a non-monotonic function. 
The requantized level density, defined as 

D T .(E,X)= J £S(E-E V {T* } X)) ) (30) 

can be rewritten in the form 

D T * (E, X) = p T * (E, X) ]T 6(Nt. (E, X)-v- 1/2) 

V 

+00 

= p T *{E,X) J2 (-l) k exp[2irikN T 4E,X)\, (31) 

with pt* (E, X) standing for the parametric level density given by the truncated trace for- 
mula. It is immediate to see that averaging Eq. (|3l| ) yields 

(Dr.{E,X)) EtX = (p T *(E,X)) E!X = p^\E,X) , (32) 

equivalent to the usual level density, easily obtained from the Weyl formula. 

The advantages of this procedure become evident once we analyze the spectral fluctua- 
tions. Let us study the 2-point spectral correlation function for Dt*(E) in order to compare 
with the results of the preceding sections. It is useful to write the density correlation function 
in the form 



K(n,X) = (d t * (£ + §,X+!) D t * (e- 2 X - f )) - (p w -') 2 • (33) 



n ^ x\ „ / „ n ^ x\\ / w„..,\2 

The averages are defined as in Section III. Thus, 



K(Q,X) = (p T ,[E + ^,X + ^)p T , ( E -^,X-^ 



x Y, (-l) kl ~ k2 expUm 

ki,k 2 \ 



n — x\ , / n — x 

/,,A> [E+-,X + -j -h 2 N T * \E--,X-- 



E,X 



-(p Wc ^) . (34) 

One can recover the standard semiclassical result discussed in the previous section by con- 
sidering the simplest term in the double sum, namely, k = ki = = 0: 
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K k=0 (tt,X) 



fi 
Pt* 



^ VL — X\ fl 



2 ' 2 



E,X 



(35) 



This is the same calculation as presented in the previous section, except that now the expo- 
nential cutoff is substituted by a sharp truncation at T*. This causes a spurious oscillatory 
term which only becomes negligible in the asymptotic regime where aX 2 /2h ^> A. In 
other words, the semiclassical result only becomes independent of the cutoff at the asymp- 
totic limit, where it coincides with RMT. In Ref. [|ITJ, this term with k\ = k 2 = was 
called "diagonal approximation". This should not be confused with the standard diagonal 
approximation discussed in the previous section. 

We are left to discuss the terms with ^ ^ and k 2 ^ 0. It is easy to see that the ones 
with ki k 2 give a negligible contribution to the correlation function. Indeed, dividing 
the spectrum for the unfolded system (|6]) into a smooth and a fluctuating part (N T * = 
jyrweyi + N fij we have ^Weyi^^x) = E / (27ih) d . For this part, the average @ furnishes 



^exp |2tu A;xiV Wcyl (e+^,X 

(2TTh) d . \(k 1 -k 2 )5E 



-)-k 2 N^ ( E --,X-- 
2 J V 2 2 



E,X 



{h - k 2 )5E 



(2irhY 



ship! - k 2 )5N] 
(h - k 2 )5N 



(36) 



which is negligible if the averaging region 5E extends over many states, unless k\ = k 2 . 
Thus, 



K k ^°{Q,X) 



°~ E Ti ex P \ 2mk (^i ~ fi 2 )p Weyl l 2 , X u X 2 ) (37) 



4vr 2 m x dSh t^o p 



(evaluated at Q\ = —Q 2 = Q and X\ = —X 2 = X), where k = k\ = k 2 and 

$ k {n 1 ,n 2 ,x 1 ,x 2 ) = (expfai^NfiiE + n^x + Xx) - n£1{e + n 2 ,x + x 2 



}) 



E,X 



(38) 



One is now left with the formidable task of averaging (38), for which there is no clear con- 
trollable strategy. Bogomolny and Keating [TjJ proposed a kind of Gaussian ansatz, which 
in practice simplifies the average by making (exp[iG(E)]) = exp[— (G 2 (E)) /2]. Although 
appealing, this ansatz is harder to justify classically than a similar passage in the previous 
section. Moreover, it is not clear how to implement the averaging in a consistent way. In 
Appendix C, we justify the validity of the Gaussian assumption (for large energy differences 
as compared to A) using solely the diagonal approximation Based on that, we have 



$ k (n 1 ,n 2 ,X 1 ,X 2 )=exp< -27r 2 k 2 { N T *(E + £li,X + X 



N T *(E + n 2 ,X + X 2 ) 



(39) 



E,X. 



There are at least two straightforward ways to evaluate the average in (^9|), both giving the 
same final result. The first one is by direct subtraction of the cumulative densities, followed 



11 



by the use of the diagonal approximation. The second one has the nice property of putting 
in evidence some of the delicate points of the semiclassical averaging procedure, and will 
be discussed now. An important point to have in mind, however, is that the average in 
the exponent of Eq. ( |3~9"D is just the usual number variance E 2 (f2) when X\ = X 2 . If we 
were interested in X = correlations at Q <C A and simply substituted into Eq. ( |39| ) the 
exact RMT expression for £ 2 (f2), we would obtain an incorrect result for the level density 
correlation function. We interpret this fact as a clear indication that the Gaussian ansatz 
of Ref.[10], which neglects correlation between orbits, is problematic for O/A <C 1. This 
limitation will appear more explicitly below. 
Let us start writing (|39|) as 



<t> k (n 1 ,n 2 ,x 1 ,x 2 



A(n 1 ,n 2 ,x u x 2 ) 

An 



with 



In A(O x , n 2 , X h X 2 ) = 4tt 2 (nV {E + Q u X + X 1 )N^ (E + Q 2 , X + X 2 



E,X 



(40) 



(41) 



and A = A(0, 0,0,0). Then, we evaluate A using the diagonal approximation (see for 
instance p3|). 



N T .(E,X) 



_9_ 

jc 2vr 2 



E,X 



, 1 

dt- 
t 



2vr 2 



T 



(42) 



The lower bound of integration is of the order of the shortest periodic orbit, beyond which 
we cannot even extrapolate the uniformity principle. Recalling that in our case the energy 
of the shell is identified with its volume, we obtain, from purely geometrical considerations, 
that 



T = |^ oc E-WV . 
dE 

If we now choose T* = Kin (t# is the Heisenberg time), Eq. ([42"D becomes 



N KtH (E,X) 



i--im 



-) 



) + constant . 



(43) 



(44) 



It is interesting that we recover the RMT result only in the limit where the phase space 
dimension d — > oo. (Incidentaly, the above formula accounts for the spectral rigidity of the 
integrable limit of d = 1.) 

Applying now the diagonal approximation to the numerator of (|40"D , we write 



\nA(tt 1 ,tt 2 ,X 1 ,X 2 ) = 2g Re f dt-exp 

Jt t 



VLo)t 



ajXi - X 2 )H 

2h 2 



which can be expressed in a compact form as 

lnA(fi 1 ,fi 2 ,X 1 ,X 2 ) = 2g Re [E^t) - E^T* 



(45) 



(46) 
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with £ = — + a(Xi — X 2 ) 2 /2h 2 and Ej denoting the exponential integral [14 



Recalling that T* is of the order of h/A and r is a much shorter time scale, we have |£T*| ^> 1 
and |£t| < 1. As a consequence, Ei(£T*) is exponentially small and Ei(£r) can be expanded 
up to first order fifM. The result is 



\nA(fl 1 ,Q 2 ,X 1 ,X 2 ) 



2g Re [-7 - ln(f r)] 
,'T' 

25- Re -7 + In 



ln(£T* 



(47) 



where 7 is the Euler constant. Collecting the above results into fl39|), we see that the first 
two terms in the r.h.s. of (f47"D are cancelled by lnA , and 



^ k (n 1 ,n 2 ,x 1 ,x 2 ) 



n 2 ) + 



a(X 1 -X 2 ) 
2h 2 



e 7 T* 



-2<?fc 2 



(48) 



Hence, the sum over k in (|37|) is an asymptotic series in inverse powers of £T*. To obtain 
the parametric generalization of Bogomolny and Keating result | I0[] , we cutoff this series in 
the first term and choose T*e 7 = tjj/2, obtaining 



m = cos(2^A) 

v ' ; 2A 2 



A 



4/0 



(49) 



depending on the symmetry class 



which, after proper rescaling, is exactly (|13D or 
labelled by (3 = 2/g. 

It thus appears that requantization has reproduced precisely the asymptotic form of the 
correlation function calculated exactly within RMT if we postulate the same factor of \[2 
for the change in X c from GOE to GUE. However, in taking a sharp cutoff for the periodic 
orbit sum, we add the spurious oscillatory term 



gT* cos ( ^1 



V h 



i>J exp ^~ 



aX 2 T* 
2ti 2 



2n 2 H 



(50) 



to (p^) with 7] = [above, = ta.n~ 1 (2hQ/aX 2 )]. This will certainly be negligible if 
aX 2 T* /2h 2 ^> 1, but we cannot push aX 2 /2h down to the averaged level spacing. Indeed, 
our deduction of ( f!9"D does not hold at this level either, since there one should not neglect 
the upper limit of (|45|). 

These conclusions are maintained if we substitute Nt* by the integral of the smoothed 
density of section [Hj in the requantization of Bogomolny and Keating. The steps in the 
calculation of the correlation function remain the same, so that now K k=0 (Q, X) is exactly 
(f28|) with a finite 77. The expression for 

/°° 1 ( r 11 

dt- exp \i{p. x - Q 2 )t/h - [a{Xx - X 2 ) 2 /2h 2 + 2r]/h\ t} 



2g Re E x 



t 

-i(fi! 



Q 2 )r/h + a{X 1 - X 2 ) 2 r/2h 2 + 2r]r/h 



(51) 



now holds, in principle, for any value of Qj or Xj, including zero. We can now use the same 
approximation for the exponential integral function as in ( f4"TD to obtain 
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$ k {n 1 ,n 2 ,x 1 ,x 2 ) 



a(X 1 



2tf 



2r)/% 



h_ 

2rj 



-2gk 2 



(52) 



which leads to (^) when we choose rj = A/2n and assume Q and aX 2 /2h to be much 
greater than the mean level spacing. In this sense, the energy smoothing with 7] corresponds 
to the cutoff time T* (r)T* ~ h). 



V. FINAL DISCUSSION AND CONCLUSIONS 

The Gutzwiller series does not converge in the limit of small smoothing. Indeed, it can 
be argued that the exponential smoothing used in Section |T| is insufficient, but Gaussian 
smoothing does work and leads to equivalent results. The application of the periodic orbit 
theory for nongeneric long-range oscillations of the energy level correlation functions was 
pushed in the last decade to a range where contact could be made with the universal regime 
of Random Matrix Theory. This involves the uniformity principle for long periodic orbits, 
as well as the diagonal approximation. In our study we have used the semiclassical trace 
formula to show the universality of the level density parametric correlation function for 
classically chaotic systems. This universality manifests itself after a proper rescaling of the 
correlation function variables. Within the range of validity of the semiclassical formulation 
and in the diagonal approximation, we have obtained the same functional dependence on 
energy and external parameter found by the method of RMT. 

From a purely utilitary point of view, it may appear unnecessary to rederive RMT results 
within a semiclassical theory. Yet, it is only in this way that we can calculate the arbitrary 
parameters in RMT, as well as derive the nongeneric long-range oscillations characteristic 
of the individual systems that we wish to measure. In the present case, we obtain the result 
that the unfolded mean square parametric velocity diminishes by a factor of \/2 when the 
time-reversal symmetry is broken. 

The semiclassical description of parametric correlations relies on the ansatz proposed 



in (^), which is based on the central limit theorem and leads to (27). To demonstrate 
rigorously the validity of (^) for a generic classical chaotic system, a systematic study of 
higher moments of Q y is required, which is a quite difficult task. Notwithstanding, there is 
solid numerical evidence p|,|2lH to support the proposed Gaussian ansatz. As a consequence 
of (p?[), aX 2 /2h will always appear in the same footing as Q in the semiclassical approach. 
In general, this is not the case in the RMT, as one can see in Section |I[ Only after 
linearizing the "effective action" (what we do to generate the asymptotic expansion), we see 
the semiclassical structure emerging. 

It seemed that one would have to proceed beyond the diagonal approximation to recover 
the oscillatory behavior which RMT predicts at short scales. However, the remarkable 
requantization scheme advanced by Bogomolny and Keating shows that it is only necessary to 
feed in the discreteness of the quantum spectrum to obtain an expression for the correlation 
functions that extrapolate to the correct oscillatory RMT result in the asymptotic limit of 
large f2 2 + a 2 X A /Ah 2 . It is also important to stress that only diagonal information about 
periodic orbits enters into this result. In other words, nothing is said about correlations 
among orbit actions, so that these may only affect higher terms in the asymptotic expansion. 
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In fact, to obtain higher-order corrections matching the exact series, one would certainly 
need to introduce information about inter-level correlations as well, which is apparently 
beyond the capability of any present semiclassical approach, specially in the case of arbitrary 
parameter variations. 

The intrinsic limitation of the semiclassical method to cover small energies at the quan- 
tum scale is remarkably manifest in the parametric correlations. The lack of accuracy in 
energy ranges Q ^$ A imposes a limitation in the accuracy of parametric correlation function 
for X < X c , as one can see from Eq. (p8|). In the standard derivation of Sec. ITTTL this is a 



direct consequence of the necessity of smoothing the level density. Even in the requantiza- 
tion scheme, although not explicitly, the same problem occurs, since both sharp and smooth 
cutoffs in the Gutzwiller series do not assure convergence to the actual eigenvalues. 

Finally, remark that in this work we did not attempt to investigate deviations from the 
universal, RMT behavior due to large scale structures in the spectrum which can be ulti- 
mately related to short periodic orbits of the system (see, for instance, the supersymmetric 
treatment used in Ref. |24]] for finite conductance disordered systems). Again, the requanti- 



zation scheme seems to be a good starting point for such systematic studies from a purely 
semiclassical point-of-view. Unfortunately, at present we only know how to proceed safely 
by restricting ourselves to the energy range where the diagonal approximation is accurate. 
There is a possibility that the resummation technique could, in principle, extract additional 
information from the A function, but this work has still to be done for generic systems. This 
is one of the major challenges of the semiclassical theory. 
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APPENDIX A: RMT ASYMPTOTIC EXPANSIONS 

Here we carry out the asymptotic expansion of k(u,x) to lowest order in x~ 2 and oj~ 2 
starting from the exact RMT expressions (f|), (f|), ([TT|) , and (0). Although, for the GUE, 
one could start from (|H1|), we decided to use the integral expression to illustrate the method. 
Essentially, we have to separate the contributions coming from the two maximum-amplitude 
points. For the GUE case, we begin by expanding F GUE around A = Ai = 1: 

1 roc roo 

k G ± E (uj,x)~-Re / ds ds\ exp[— (n 2 x 2 — muj + ){s + Si)] 

2 jo Jo 

(Al) 



-Re 
2 



(tt 2 x 2 — muj + y 



where s = 1 — A and Si = Ai — 1. Notice that the integrand in ( |A1| ) is only appreciable for 
s, Si ~ 0(x~ 2 ) and, since we only need the lowest order in x~ 2 , we can neglect quadratic 
terms in s. Equation (|A1|) , however, is still correct even if x — 0, provided that keep lu 1. 
Next, we expand F GUE around A = — Ai = — 1: 
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1 roo roo 

k^Yfiu, x) « -Re e 2t7TUJ / dr / ds t exp[-vrV(r + s x ) - inco + (r - s x )\ 

2 Jo Jo 



1 / e 2inuJ \ 
« »R* | 2 2 • + , 2 , (A2) 

where r = A + 1. This result is also valid for x = and to ^> 1. Notice moreover that, 
contrary to the previous contribution, this one is oscillating in uj on a scale 0(1) and therefore 
contains information about the discreteness of the spectrum. We can add ( JA1| ) to ( |A2p to 
obtain Eq. (0). 

The GOE case is similar, but slightly more complicated due to the additional structure 
in the integrand. Beginning with A = Ai = A 2 = 1, we can expand F and the rest of the 
integrand around this point to get 

r°° roo (2s)(s + Si + S 2 f 



roo roo roo 

fcfff (w,x) « Re / ds dst ds 
Jo Jo Jo 



' 2 (2ssi + 2ss 2 - 2 Sl s 2 + sj + si + s 2 ) 2 
x exp[-(vrV/2 - mcj + )(s + si + s 2 )] • (A3) 

where Sj = Aj — 1. This integral can be evaluated by using spherical coordinates in the first 
quadrant, giving 

,GOE 



k%X?(w,x) « Re 



(7T 2 X 2 /2 — i-KU+y 



r /2 , , T /2 ,„ sin 20 

x / / Gt0 



Re 



o Jo [1 + sin2#(cos</> + sin0) + (cos 20 + 1) sin cos 0] 2 

(A4) 



(ti 2 x 2 /2 - mu+) 2 



(Here we have resorted to a numerical method to solve the double integral.) Finally, we can 
expand around the other maximal point, which is A = — Aj = — A 2 = —1, and obtain 

(2r)(-2) 2 



k^fAu^x) « Re |e W [°°dr [°° d Sl f°° ds 2 ^ 
[ Jo Jo Jo (— 



_ 4 )2 

x exp[-(7r 2 x 2 /2)(r + si + s 2 ) - muj + {r — S\ — s 2 )} j 

1 / 21TTU1+ \ 

« - Re 7^77 Z77 • ( A5 ) 

Since this contribution is of higher order than we should only keep the oscillating part 



and neglect the rest. With this in mind, we arrive at Eq. (14). 



APPENDIX B: PARAMETRIC VARIATION OF PERIODIC ORBIT 

It may be surprising that we can integrate the parametric velocity Q 7 in (0), obtained 
from classical perturbation theory to yield the exact result: 
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S(E,X 2 )-S(E,X 1 ) = J x *dX jf ' dt — (t{t),X) (Bl) 

where £ = (p, q) are vectors in the classical phase space. Thus there is no limit on the period 
of the orbit if we keep the right order of integration. We can obtain the expression above by 
embedding the one parameter family of Hamiltonians H(£,X) into a single Hamiltonian 7i 
in a phase space that is expanded by two more coordinates. Hence, we add the parameter 
X itself and a conjugate variable Y, defining the Hamiltonian so that 

H^,X,Y)=H^,X) (B2) 

at each point. Then, Hamilton's equations determine that X is a constant of the motion, 
whereas the equations for £ are unaltered, so that the energy of the original systems is still 
constant, equal to E. However, the periodic orbits of the original systems now correspond 
to helicoidal orbits such that 

T T /~) TT 

AY = Q,(E,X) = J Q J dtYm,X)= j o 7 dt — (at),x) . (B3) 

Conservation of E and X implies that each of these helicoidal orbits lies within a two- 
parameter family in the extended phase space. Fixing the energy E, we thus obtain a two 
dimensional surface, along which 



j> dq 



p + j>dXY-j>dtH = (B4) 
for any reducible circuit (as a consequence of the Poincare-Cartan theorem [25]). The last 



integral cancels, because the full energy has been chosen identical to the energies of the 
original system, a constant along the surface. Picking the circuit so as to connect two 
helicoids (corresponding to periodic orbits) with different parameters, we obtain 



I dq • p - I dq • p = f ' dX Y , (B5) 

recovering (|HT|). 

It is important to note that the action difference refers to a continuous family of periodic 
orbits. In the case where there is a set of multiply symmetry-connected orbits, (Bl) can 
also be used if the change of parameter does not break the symmetry. 



We also point out that (|B1| ) refers strictly to an action difference at constant energy. Of 
course, it is possible to choose different one-parameter families of periodic orbits within the 
two-parameter family in our problem. Goldberg et al. [[J] choose constant volume for the 
shell instead. This coincides with constant energy for unfolded dynamical systems. 



APPENDIX C: JUSTIFICATION FOR THE GAUSSIAN ANSATZ 



In this appendix we motivate the Gaussian ansatz used in Ref. [TO] using arguments based 
solely on the diagonal approximation. Let us start by Taylor expanding the exponential term 
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of of Eq. fl38|), thus recasting the original exponential averaging problem into an averaging 
of increasing powers of (5N^[) n , with 5Ng, = N{L(E,X) - N#.(E',X'). 

The lowest-order term of the Taylor expansion (n = 1), taken within the validity range 
of classical perturbation theory, can be written as 



N^i (E + n 1 ,x + x 1 )- n{!. {e + n 2 ,x + x 2j 



irSJE + Q u X + X X ) 



E,X 




exp 

which clearly vanishes, since 

JrSj(E,X)/h 



Tl 



exp 



irS 7 (E + n 2 ,X + X 2 ) 



(Cl) 



e irT 7 (_B,X)ni/^+irQ 7 (£;,X)Xi/^ _ e irT 1 {E,X)il 2 /h+irQ 1 {E,X)X 2 /n 



E,X 



E,X 







(C2) 



due to the very rapid oscillations of the term in S 1 /h in the averaging interval 5E. 

Generally, in order to energy average the n-th power of 8N^{ we have to calculate terms 

like 



An) 
>,7> 



exp 



^{nS^ + r 2 S r2 + 



+ r n S r 



(C3) 



In fact, neglecting repetitions (see comment bellow), we conclude that w 0, whenever n is 
odd, as in (Cl). Moreover, for even values of n, we can always find a set of actions that cancel 
each other, yielding I/" 7 \ ^ 0. Due to phase-space restrictions, the largest contribution 

will come from sets where the actions are grouped pairwise. This allows ((6N^) n ) to be 
factorized into powers of ((SNj,[) 2 ), leading to the Gaussian formula used in Ref. [K|. Note, 
however, that our argument relies on 



-(riS 1 ^ + r 2 S l2 H V r n S ln ) > 1 



(C4) 



for any distinct combination of trajectories. This condition is not necessarily satisfied when 
n > 1. 

Finally, in a first inspection, the inclusion of repetitions would seem to spoil our argu- 
ments in favour of the Gaussian ansatz. However, using again phase-space considerations, 
it is simple to see that the number of terms obtained by cancelling actions using repetitions 
is much smaller than the number of simple pairwise cancellations of primitive orbits. Thus 
the effect of the repetitions is not expected to be important. 
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